arXiv:1508.07743v3 [math-ph] 19 Feb 2017 


On Poincare’s generating function and the 
mid-point rule 

H. Jimenez-Perez, J.-P. Vilotte, B. Romanowicz 
February 21, 2017 


Abstract 

The use of Liouvillian forms to obtain symplectic maps for construct¬ 
ing numerical integrators is a natural alternative to the method of gen¬ 
erating functions, and provides a deeper understanding of the geometry 
of this procedure. Using Liouvillian forms we study the generating func¬ 
tion introduced by Poincare (1899) and its associated symplectic map. 
We show that in this framework, Poincare’s generating function does not 
correspond to the symplectic mid-point rule, but to the identity map. 
We give an interpretation of this result based on the original framework 
constructed by Poincare. 

Additionally, we construct families of Liouvillian forms generating, 
as a particular case, the mid-point rule, showing that their structure is 
different from the differential of the generating function used by Poincare 
for studying periodic orbits. 


1 Introduction 

A symplectic integrator for a Hamiltonian system is a numerical method which 
preserves the structure of the Hamiltonian vector field. From the available tools 
for constructing symplectic maps, the method of generating functions has been 
a cornerstone to understand the links between the geometry and topology of the 
phase space of Hamiltonian mechanical systems. Its relevance resides in the fact 
that this method is used for theoretical and applied results in contrast to other 
methods which are only valid for either theoretical or numerical applications, 
but not both. 

In his famous Les methodes nouvelles de la mecanique celeste [II], Poincare 
develops the theory of integral invariants with applications to the study of pe¬ 
riodic orbits in celestial mechanics. Poincare constructed an exact differential 
1-form defined on closed orbits with prescribed fixed period T > 0, such that 
its exterior differential gives the canonical symplectic form on the phase space. 
Since the orbit was periodic he considered a section (the Poincare’s section) and 
a non-trivial map defined on the section (the Poincare’s map) such that the pe¬ 
riodic orbit corresponds to a fixed point of the map. The imposed condition for 
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the fixed point was that the first-return map must be non-reversing. All these 
hypotheses will be important for the interpretation of the results described in 
this paper. 

From the numerical point of view, generating functions are used to con¬ 
struct numerical algorithms preserving the main geometrical properties of the 
phase space. These properties come from a geometrical structure known as the 
symplectic form, naming the numerical algorithms as symplectic integrators. 
The term symplectic was coined by Weyl in 1946, however, some geometrical 
properties were already used by Euler, Lagrange, Jacobi, Hamilton, Poincare, 
among many other mathematicians and physicists of the last three centuries. 
Symplectic integrators arose after the construction of the main computers, and 
what seems to be the first documented reference, is the paper of De Voge- 
leare [13] in 1956. In the early 80’s several works from Ruth [12], Channell [1], 
Menyuk [10] and Kang [2] on symplectic integrators aroused the interest of the 
numerical community because of their excellent behaviour and stability for long 
time computations. In the second half of the 80’s, symplectic integrators using 
generating functions were systematically studied by Feng Kang and co-workers 
[2, 8, 3, 4]. In all those papers, the Poincare’s generating function is systemat¬ 
ically associated to the symplectic mid-point integrator. However, there is no 
formal demonstration or construction of this correspondence. 

Recently, a more general method for constructing symplectic maps has been 
developed using Liouvillian forms [5]. This method gives a better perspective 
of the role played by Lagrangian submanifolds in the search of symplectic maps 
for constructing symplectic integrators. Moreover, it has been proven that Li¬ 
ouvillian forms coupled with Hamiltonian vector fields produce a function which 
gives the variation of the Hamiltonian function along the Liouville vector field 
(the symplectic dual of the Liouvillian form) [6]. This function has an intrin¬ 
sic symmetry, since it also gives information on the variation of the Liouvillian 
form along the Hamiltonian vector field. Moreover, the difference between this 
function and the Hamiltonian function is a sort of “misfit” function called the 
elementary action, whose integral with respect to time is the action integral of 
classical mechanics [9]. These results yield another point of view on symplec¬ 
tic integration since we will be interested in the search of the Liouvillian form 
which produces the function as close as possible to a constant function on the 
integral lines of the Hamiltonian vector field. The direct consequence is that the 
elementary action will be constant on the solutions of the Hamiltonian vector 
field, satisfying automatically Hamilton’s principle of minimal action. Since we 
work on autonomous Hamiltonian systems, the elementary action must be a 
function of the Hamiltonian but not of time, which is the common symmetry 
(hrst integral) to all autonomous Hamiltonian systems. 

In order to produce symplectic integrators, a study of the structure of Li¬ 
ouvillian forms with constant coefficients has been achieved in [5]. The main 
result is that all the one step implicit symplectic integrators are given by the 
mid-point rule plus a Hamiltonian component which vanishes in the limit h —>■ 0, 
where h is the time step. Consequently, it is necessary to study the Poincare’s 
form, i.e. the differential of the so called Poincare’s generating function, for 
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comparing its structure with the differential form related to the mid-point rule. 
We develop this analysis by applying the method of Liouvillian forms [5] to 
Poincare’s differential form, with a set of remarks in order to clarify some mis¬ 
understandings about both objects. In addition, we find a non-trivial family 
of Liouvillian forms (on the product phase space) associated to the mid-point 
rule. In contrast to the claim found in the papers of Kang and his collaborators, 
our results associate the differential of the Poincare’s generating function to the 
identity map. This is not a surprise since both mid-point rule and Poincare’s 
differential form belong to different families of minimizers of the action integral. 
Indeed, the mid-point rule minimizes the action along a path with different fixed 
boundary points. Meanwhile Poincare’s differential form minimizes the action 
integral along a periodic closed path with prescribed fixed period T > 0, char¬ 
acterized on a Poincare’s section by a fixed point. We interpret this result using 
the original hypotheses stated by Poincare in [11]. 

2 The method of Liouvillian forms 

In this section we give some general facts about the method of Liouvillian forms 
and generating functions, and we briefly recall how the procedure is applied to 
construct the symplectic integrators. 

We consider the momentum phase space of a mechanical system as a sym¬ 
plectic manifold (M = T*Q,w) with the standard symplectic form satisfying 
wo(-, •) = (•, Jq-)- Here (•, •) is the standard inner product and Jq is the canoni¬ 
cal complex structure or canonical symplectic matrix given by 

Jo = ( o” ) , In,Or, G M„xn(M). (1) 

A Liouvillian form on a symplectic manifold (M, uj) is a 1-form 9 on M such 
that d9 = Cl). In this case we say that {M,d6) is an exact symplectic manifold. 
If 6 is the tautological 1-form we call it the Liouville form, in other cases we 
call it a Liouvillian form. A Lagrangian submanifold A C M in a symplectic 
manifold (M, w) is a submanifold where the symplectic form vanish identically 
w|a = 0. 

Lagrangian submanifolds are the integral submanifolds of the Liouvillian 
forms where the integrability is defined in the following way: given a Liouvillian 
form 9, the Lagrangian submanifold A is an integral submanifold of 9 if 0 |a = 0 
and Cl) IA = d9\A = 0. Along the Lagrangian submanifold A it is possible to 
decompose the tangent space of the symplectic manifold in a direct sum of 
vector subspaces T^M = t4A 0 T^A for every z G A. I 4 A is the vertical space 
of A at z in T^M and we can recover it using the complex structure J by 
t4A = JTzA. It means that for every z G A and v G T^A the vector Jv G I 4 A. 

Geometrically, a Hamiltonian mechanical system corresponds to a Hamilto¬ 
nian vector field Xh on the symplectic manifold (M, w) representing the mo¬ 
mentum phase space. The dynamical system is given by the equations of motion, 
also known as Hamilton’s equations z = X//(z), where Xh = JqVH, V being 
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the canonical gradient operator. A first order (generic) approximation is given 
by the Euler map 

z.h = zo +hXniz), z,Zh,zoGM, h G R+, (2) 

where z is an element to be determined. In order to find z G M, such that 
(2) becomes a symplectic map, we use the classical framework for constructing 
symplectic maps using generating functions. 

Consider the product manifold M = Mi x M 2 of two copies of the phase space 
with the bilinear form wq = ttiOJi — 7rja;2 induced by the canonical projections 
TTi : M —>■ Mi, i= 1,2. The couple {M,ujq) becomes a symplectic manifold of 
dimension 4n [9]. 

Given a Euclidean structure (•,-) 4 n on T^M for m G M, we associate the 
complex structure J to ujq(-, ■) = (•, J-), where 

J=( Jo 02n \ g M2„x2n(M) 

\ U2n —Jo ) 

and Jo as given in (1). 

In this framework, the method of Liouvillian forms consists in finding a 
couple of related Liouvillian forms Bi on the manifolds {Mi,u!i), i = 1,2, such 
that the integral submanifold A of the induced form 9 = — 7r202 on (M, Wq), 

is adapted for constructing a numerical scheme. The method follows four steps: 

1. Eind a suitable Liouvillian form 0 G T*M constructed by the pull-back of 
two Liouvillian forms Bi on [Mi,uji), i = 1, 2 in this way B = tt^Bi — 71202- 

2. Obtain the expression of the vertical bundle EA to the Lagrangian sub¬ 
manifold A C M associated with B (writing its coefficients in local coor¬ 
dinates) . 

3. Project the tangent space of the Lagrangian submanifold TA C TM into 
the tangent space TM oi the original symplectic manifold {M,lo). We 
find the local expression of TA using the equation 0 = 0 and the complex 
structure TA = J^VA. We consider the sum T7ri(TA) -|- T7r2(TA) as 
vector subspaces of TM and we define the implicit map 

p{zo,Zh) = T7ri(T(^^,^,^)A) -k T7r2(T(^„^^,^)A) 

4. Verify that we recover the tangent space TM when we evaluate TA C TM 
on the diagonal. Eormally we check if T^M = T7ri(T(2; z)A)0T7r2(T(2; z)A) 
is satisfied. This implies that p becomes the identity map p(z, z) = z. In 
this case the generalized implicit Euler map 


Zh — zq -t- hV// o p(zo, z/j). 


(3) 


is symplectic. 
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The interested readers are referred to [5] for the details of the construction. 

We can check that the Liouvillian form 0 on {M,u>q) is adapted for the 
construction of a symplectic integrator if the induced implicit map becomes 

p(zo,z/,) = i(zo + z/i) + b(z/i - zo), (4) 

where b £ M2rax2n(R) is a Hamiltonian matrix [7, 5]. In particular, the case 
b = 02n corresponds to the mid-point rule. In the next section we will apply 
this construction to the generating function introduced by Poincare in [II]. 


3 The Poincare’s Generating Function 

In the third volume of Les methodes nouvelles de la mecanique celeste [11], 
Poincare introduced the exact 1-form 

^ E {(Q - + p) - (P - P)d{Q + q)} (5) 

looking for periodic orbits bifurcarting from a prescribed periodic orbit of period 
r > 0. In expression (5) variables (q, p) are positions and conjugate momenta 
in the phase space at time t and (Q, P) are positions and conjugate momenta at 
time t + T. We denote them by Zq = (q, p) and Zh = (Q, P). The corresponding 
function S' : A C M —> R is called the Poincare’s generating function. This 
exact form was rediscovered by Feng Kang and his collaborators when they were 
studying the construction of symplectic integrators using generating functions. 
Kang’s group interprets this 1-form as the linear mapping 

pdq -f PdQ ^ [(Q - q)d(P -f p) - (P - p)d(Q -f q)] (6) 

given by the matrix 

Jo,l2n s M2nx2ri(R), (7) 

where Jq is the canonical complex structure (1) and l 2 n is the identity matrix 
of size 2n. Given a generating function u : A —>■ R on a Lagrangian submanifold 
Ac M with local coordinates w = w{zQ,Zh), they systematically associate the 
numerical method 



du f Zh+zo 

Zh — Zq = -Jq-^ 

ow 


( 8 ) 


to the exact form (5) [2, 8, 3, 4], and consequently the symplectic mid-point 
rule with the Poincare’s generating function. However, the definition of the 
Lagrangian submanifold given by the 1-form and the Hamilton-Jacobi equation 
were not formally established. Note that the map (8) is a time-1 map. 
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The mapping given by the matrix a in (7) is not unique. In particular we 
can consider the following matrices 
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On 

It 

2^n 

On 

1 / 

2-^n 



/ In 

On 

In 

On 

\ 



i/n 

On 

^In 

On 


1 

On 

In 

On 

In 


ai = 


-In 

On 

In 

On 



On 

In 

On 

-In 



V 

On 

-In 

On 

In 



V -In 

On 

In 

On 

/ 


generating the same mapping for the differential 1-form ( 6 ). They have a very 
different structure. 

The matrix a given by Kang is motivated by the use of the identity matrix 
l 2 n and the complex structure Jq acting in a natural way on vectors of the 
product manifold. In their interpretation, the left hand side of ( 8 ) corresponds 
to the first “row” of the matrix a, and the argument in the right hand side, 
corresponds to the second “row”. However, such an interpretation loses meaning 
if we write the same mapping using matrices ai or 04 . 


Remark 1 The fact that we obtain the same Liouvillian form for very different 
linear transformations on M comes from a symmetry induced by the construc¬ 
tion of the product manifold {M,loq) which do not correspond to symmetries on 
the original manifold {M,u}). 


Given a Liouvillian form 6 on the product manifold M, we can apply the 
method above defined to construct an implicit symplectic integrator. Let us do 
that to the 1-form (5). 

Step 1) This step does not apply since we already have a Liouvillian form. In 
fact, our goal is to prove that Poincare’s 1-form is not suitable for constructing 
a symplectic integrator. For convenience we rewrite the Poincare 1-form (5) as 


e = i[(p-P)dq+(Q-q)dp + (p-P)dQ+(Q-q)dP]. (9) 


Step 2) The local expression of the vertical space I 4 A at z = (q, p, Q, P) G A 
to the Lagrangian submanifold is determined by the coefficients of the Poincare’s 
1-form 9. They correspond to the algebraic equations 


q = |(p-P), p = |(Q-q) 

Q = i(p-P), P = i(Q-q). 


( 10 ) 


where we consider q, p, Q, P : TM —>■ R as real functions on the product mani¬ 
fold. 

Step 3) The local equations of the tangent space TzA are obtained by the mul¬ 
tiplication of the complex structure with the vertical space T^A = J^14A. Using 
the expression (10) the complex structure J maps (q, p, Q, P) >->• (—p, q, P, —Q) 
whose induced implicit map is given by 


p(q, p, Q, P) = (-P, q) + (P, -Q) = (P - p, q - Q)- 


6 







Figure 1: The original framework of Poincare. The direction of the perturbation 
is perpendicular to the flow’s direction. 


Substituting (10) yields unexpectedly: 

p(q,p,Q,P) = (0,0). (11) 

We consider that Xif(0,0) = 0 as Poincare did. 

Step 4) Composing the generalized Euler integrator with the map p we obtain 

= zo + hX//(0, 0) = zo (12) 

which corresponds to the identity map. 

This fact is consistent with Poincare’s hypotheses, since he constructed the 
generating function to be defined on periodic orbits of prescribed period T > 0 
such that if zg = (q, p), '^h = (Q,P) then z^ = zg was a fixed point on 
a Poincare’s section but such that the orbit was not trivial (see Figure 1). 
Moreover, this hypothesis considers that for T variable, limT->g dS = 0 and 
the function S goes to a constant S ^ So which, for simplicity, he considered 

5'g = 0. 

We have proven the following: 

Proposition 3.1 The symplectic map associated to the Poincare’s form (5) 
under the method of Liouvillian forms corresponds to the null map (11). If 
X^f(0,0) = 0 the generalized Euler scheme corresponds to the identity map. 

Remark 2 Another way to see that Poincare’s form (5) corresponds to the iden¬ 
tity map is by checking the kernel of 6 as a linear form. This directly gives the 
equations (Q - q) = 0 and (P — p) = 0, producing the identity map Q = q and 

P = P 

3.1 A path of Liouvillian forms 

To characterize the relationship of the Liouvillian forms producing symplectic 
maps with the Poincare’s 1-form (5), we construct a path of forms connecting 
the symplectic Euler methods A and B with this 1-form. 
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Lemma 3.2 Let In, On S M„xn(R) be the identity and the zero square matrices 
of size nxn. Every element in the 1-parameter family of linear transformations 
R(j, € Diff(M,a;e) given by 

/ cos (fin 0„ sin (fin 

D COS (fin On 

On -Sin (fin On 

y -sin (fin 0 „ cos (fin 

is a symplectic rotation on T^M for every x S M. 

Proof. It is a straightforward computation checking that = I An and 

R^JR(f = J, which implies that it is an orthogonal and symplectic mapping, 
i.e. a symplectic rotation. □ 


On 

- sin (fin 
cos (fin 
On 


,0e[O,27r] (13) 


Denote hy Ei : M T*{Qi x Q 2 ) the canonical symplectic diffeomorphism 
between the product manifold {M,uj) and the cotangent bundle T*{Qi x 22 )- 
Given local coordinates (q, p, Q,P) G M and (x, X, y,Y) G T*(Qi x 22 ), the 
symplectomorphism Ei whose pull-back maps (ydx -f YdX) 1 —>■ (pdq — PdQ) 
is given by 


/ In 

On 

On 

On 

\ 

On 

On 

-In 

On 

On 

In 

On 

On 


\ o„ 

On 

On 

In 

/ 


Define the diffeomorphism : M —>■ r*(2i x 22 ) as the composition of R^ 
given in (13) and Ei. This composition determines a curve of symplectic 
diffeomorphisms and consequently, a 1-parameter family of Liouvillian forms 
= ^^^QixQ 2 on M. 

Lemma 3.3 The 1-parameter family of Liouvillian forms 9^ = 'I' 0 ^gixQ 2 given 
by 


6rj, = (cos — sin (fq)d{cos (fP -\- sin c/jp) 

— (sin(/)P — cos(fp)d{sin(fCl -\- cos^q). (14) 

connects Poincare’s 1-form, to those associated with the symplectic Euler inte¬ 
grators A and B. 

Proof. It is a family of Liouvillian forms by construction since d9^ = wq for 
every (f G [0, 27r]. To prove that it contains the specific forms, it is enough to 
compute 00 for the values (f G {0,7r/4,7r/2} obtaining 

00 = pdq + QdP, 0,r/2 = —qclp — PdQ, (15) 

^ ((Q - q)'^(P + p) - (P - p)c^(Q + q)) (i6) 
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which corresponds to the forms associated with Euler integrators A and B [5] 
and the Poincare’s 1-form, respectively. □ 


Let’s apply the method of Liouvillian forms to the family (14). We do not 
perform Step 1) since we already have the Liouvillian forms of interest. 

Step 2) The vertical space 14 A at z = (q, p, Q,P) G A to the Lagrangian 
submanifold is given by the algebraic equations 

q = cos^ (/ip — cos^sin(/)P, p = cos sin — sin^ 0q 

Q = cos0sin(/)p — sin^ ^P, P = cos^ — cos ^ sin </)q 

Step 3) The tangent space PzA = J^14A using the equation of the vertical 
space is J’^(q, p, Q, P) = (—P, q, P, —Q) with the induced implicit map 

p(q,p,Q,P) = + 

where = — sin(/)(cos^ — sine/)) and g((/>) = cos^(cos(/> — sin(/)). 

Step 4) Composing the generalized Euler integrator with the map p we obtain 
the method 


Zh = zo + hXn O p(q, p, Q, P) 


(17) 


where p is defined in Step 3). 

The generalized Euler map is a symplectic integrator if the argument of the 
discrete Hamiltonian field has the form p(zo,z/i) = i(zo -I- Zh) + b(z?i — Zq) 
where b G M 2 nx 2 n(R) is a Hamiltonian matrix [5]. Denoting Zq = (q, p) and 
Z/t = (Q,P) it is easy to prove a weak condition; if the implicit Euler map 
is symplectic the matrices of partial derivatives satisfy 9zoP + d^^p = l 2 n and 
dzhP ~ dzaP = 2b. These matrices are given explicitely by 

dp{zQ,Zh) f /(</>)/„ 0„ \ dp{zQ,Zh) f g{4>)In 0„ \ 

9zo 'v 0„ dzh \ 0„ /(</.)/„ j ■ 

The final result depends on the functions 

f{4‘)3-g{4’) = 1 — 2 sin (j} cos (^ = 1 — sin 2(j), 

f{4‘)—g{4') = — (cos^ 0 — sin^ 0) = — cos2(/). 

Both conditions are satished in [0,7r/2], if and only if ^ = 0 or ^ = 7r/2, 
where the values of b are 


b = 


1 

2 


In 



and b = i 


In On 
On In 


respectively. They corresponds to the matrices of the symplectic Euler schemes 
B and A. We have proven the following: 


9 





Figure 2: The family Off, and its projection on the subspace {Qi x Q 2 ). The set 
of Liouvillian forms whose projection gives symplectic integrators, reproduces 
the original positions and conjugated momenta on the diagonal (q = Q, p = P), 
which is not the case for the Poincare’s 1-form. 


Theorem 3.4 The family of Liouvillian forms (I 4 ) renders the generalized im¬ 
plicit Euler method symplectic if and only if tf = 0 or (j) = 7r/2, equivalently if 
and only if it is one of the symplectic Euler schemes: A or B. 

We constructed a family of implicit maps which intersects the space of sym¬ 
plectic maps in two points. The projection of this family in the configuration 
space {q, Q) is given in Figure 2. 


3.2 The Poincare’s 1-form in matricial notation 


Before presenting a family of Liouvillian forms producing the mid-point rules, 
we want to give the matricial expression of the Poincare’s 1-form in order to un¬ 
derstand the structural difference with Liouvillian forms generating symplectic 
integrators. 

Expand the family (14) and denote by Z = (q, p, Q, P)^ the coordinates of 
points in the product manifold Z S M. We can rewrite (14) in matrix form 
c?Z^[A]Z by means of the matrix 


[A] 


/ 0 
- sin^ (j)In 
0„ 

\ — cos (j) sin 4>In 


COS^ (fin 

On 

COS (j) sin (fin 
On 


On 

COS (f sin (fin 
On 

COS^ (fin 


— COS (f sin (fin 
On 

- sin^ (fin 
On 


\ 


Using the identities sin2(() = 2cossin(() and cos^ (f = (1 — cos2(f)/2 we 
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rewrite [A] as 


[A] — ij + 



/ 

0 

In 

On 

On 



cos2(() 


In 

On 

On 

On 


sin2(() / 0 -Jo \ 

2 


On 

On 

On 

In 


2 V "^0 On ; 


V 

On 

On 

In 

On 

) 



Which lets us compare the structure of Poincare’s 1-form with the Liouvillian 
forms generating symplectic integrators: if sin 20 ^ 0 , the last term on the right 
hand side is symmetric which implies that dZ^[A]Z is a Liouvillian form on 
(MjWq); but the submatrices are antisymmetric and this form cannot be the 
pull-back of a Liouvillian form on some symplectic submanifold (IV, dr]). 

Moreover, the Poincare’s 1-form corresponds to the value (j) = tt/A with 
associated matrix 


[A] = i 


Jo —Jo 
Jo —Jo 


4 An alternative family of Liouvillian forms 


In an analogous way to the 3n-parameter family introduced in [5], we intro¬ 
duce parameters (ai,(3i, ^i) i = 1,... ,n, generalizing the coefficients cos^ (p 
cos0sine/) in (14). We rearrange the elements of Otj, in order for (3i and to 
become the coefficients of elements which vanish under the differential, in the 
following way 


— (2 (PiJQi QidPi) (2 (Pidpi -f PidQj)) (f^) 

Pl^i(PiJQi P QiJpi) 'yi(QiJPi “b Pi^Qi)’ 

With this modification, the matrices of partial derivatives become: 


dp(zo,Zh) ( 

a 

Oil 


On 

1 

0 

SI 


0 


p Ui — Pi) In 

dp(zo,Zh) / 


P Oii — 

P^)In 

On 

dzh V 


On 

a 

1 

Q 

1 

3 " 

Conditions d^gP P dzhP = 

hn 

and 

;hP-dzoP = 

= 2b, for b Ha 


satisfied if and only if /3i -|- 7 ^ = 0 . 

This reduces to a 2n-parameter family of Liouvillian forms generating valid 
symplectic maps. For each pair of parameters {ai, Pi} we define a new parameter 
Si = ai — Pi for Si € [0,1]. The matrices of partial derivatives 9zop(zo, Zh) and 
5zhP(zo,z/t) becomes 


dzoP — 


— f (2 Si) In 


0„ 

(i -b S, 


)ln 


and dz^p = 


b Si) In On 

On (2 ^0 


All these maps interleave Hamiltonian maps and consequently, they induce im¬ 
plicit symplectic integrators. The particular case Si = 0 corresponds to the 
symplectic mid-point rule. We have the following: 
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Theorem 4.1 With the previous notation and setting on = Pi where i = 1, • • • , n, 
the n-parameter family of forms 


= (5 + A) {PidQi + QidPi) - (5 - Pi) {qidpi + PidQi) (19) 

+Pi{p^dQi + Qidpi + q^dPi + Pidpi). 


generates the mid-point rule under the method of Liouvillian forms. 

Proof. Let us apply the method of Liouvillian forms. The first step is to ob¬ 
tain the Liouvillian form candidates, or to verify that the proposed 1-form is 
Liouvillian. This step is fulfilled by the elements of the family (19) by construc¬ 
tion. The second step concerns the equations of the vertical spaces VAp^ to the 
Lagrangian submanifold Ap., giving 

Qi = {^ + Pi) Pi + PiPi, Pi = - - Pi) Qz + PzQi, 

Qi = — — Pi) Pi + PiPi, Pi = + Pi) Qi + PiQi- 

The third step is to recover the tangent space TAp. = JVAp^ and the construc¬ 
tion of the implicit map by the projection of TA^. leading to 

/o(q,p,Q,P) = (i(q +Q), i(p + P)), 

which is independent of Pi. 

The fourth step is the composition of the generalized Euler scheme with the 
implicit map p. We use the notation zq = (q, p) and Zh = (Q, P) for writing 

z,, = zo -b hXn (^(zo + z/i)) . 

This is the symplectic mid-point rule as we want to prove. □ 

In particular, the case /3i = 0, i = 1,..., n, gives the standard Liouvillian 
form which we can associate canonically to the mid-point rule 

do = 5 (pdq — qdp — PdQ -f QdP). (20) 

In contrast to the canonical case (20), the family of Liouvillian forms (19) does 
not accepts a description in the form 9 = TrjI'Si — 7 r 2 02 . Figure 3, shows the lines 
Si = constant for Si € [0,1] in the space of parameters (oi, /3i, 7 ^). In this space, 
the whole plane Pi "fi = 0 generates valid symplectic integrators. 

5 Conclusions 

In this paper we have applied the method of Liouvillian forms [5] for studying 
the Poincare’s 1-form^ (5) introduced in [II] and its relation with the mid-point 
rule. We showed that the classical association between these two objects is not 
the right one. This comes from the fact that Poincare’s I-form and mid-point 
rule are techniques applied to two different types of variational problems: 

^Also known as Poincare’s generating function 


12 




Figure 3: Lines Si = Ui — /3i on the plane /3i + 7 i = 0. Each line corresponds to 
one implicit symplectic integrator under the method of Liouvillian forms. 


• Poincare’s 1-form was designed for dealing with periodic orbits with pre¬ 
scribed period T > 0, it means, non-trivial loops or cycles (no boundary); 

• the mid-point rule is the simplest approximation for problems with fixed 
values at the boundary (initial and final fixed points). 

The former deals with the topology, i.e. global geometry, of the phase space 
and the latter with its local geometry. This fact becomes evident when we 
consider the limiting case producing the identity map: the loop used with the 
Poincare’s 1-form remains a loop with period T > 0 meanwhile the open path 
of the mid-point rule converges to a single point. 

We showed that the structure of Poincare’s 1-form differs drastically from 
the structure of Liouvillian forms generating the mid-point rule, and in general 
to those generating symplectic integrators. In order to better understand this 
discrepancy, we constructed two families of 1-forms. The first one is a one- 
parameter family (a path) joinning the symplectic Euler maps A and B with 
Poincare’s 1-form; this path is constructed by a symplectic rotation. The only 
elements in this family which generate symplectic integrators are the boundary 
points of the path corresponding to the Euler maps A and B. The second family 
is a 2n-parameter family of 1-forms producing symplectic maps from where any 
element of a n-parameter sub-family generates the mid-point rule. 
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